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ABSTRACT 

We present a set of numerical simulations addressing the effects of magnetic field strength and ori- 
entation on the flow-driven formation of molecular clouds. Fields perpendicular to the flows sweeping 
up the cloud can efficiently prevent the formation of massive clouds but permit the build-up of cold, 
diffuse filaments. Fields aligned with the flows lead to substantial clouds, whose degree of fragmen- 
tation and turbulence strongly depends on the background field strength. Adding a random field 
component leads to a "selection effect" for molecular cloud formation: high column densities are only 
reached at locations where the field component perpendicular to the flows is vanishing. Searching for 
signatures of colliding flows should focus on the diffuse, warm gas, since the cold gas phase making up 
the cloud will have lost the information about the original flow direction because the magnetic fields 
redistribute the kinetic energy of the inflows. 

Subject headings: MHD — instabilities — turbulence — methodsmumerical — stars:formation — 
ISMxlouds 



I. MOTIVATION 

Evidence is accumulating that star forma- 
tion follows rapidly upon mole c ular cloud 

formation (e.g. iHartmann et al.l 1200 ll and 

iBallesteros-Paredes fc Hartmanij |2007| for the so- 
lar neighborhoo d; lEngargiol a et al n2003l for M33; 
lElmegreenl l2007l in the context of M51). This rapid 
onset suggests that the clouds need to acquire high, 
non-linear density enhancements during their formation, 
since massive, finite clouds are highly susceptible to 
global gravitational collapse which could overwhelm 
small-scale fragmentation necessary for (local) star 
formation (Burkert & Hartmann 2004). Thus to under- 
stand the initial conditions for star formation, we need 
t o understand the formation of th e par ental clouds. 

iBallesteros-Paredes et al.l ()1999[ ) and IHartmann et"aLl 
()2001t ) proposed that the build-up of clouds in large- 
scale, converging fiows of diffuse atomic gas could ex- 
plain the crossing time problem, i.e. the observa- 
tion that the stellar age spreads in a large number 
of local star forming regions ar e substantially smaller 
than the l ateral crossing time ('Ha rtmann et al.l 120011 : 
IBallesteros -Paredes & Hartmann 2003). In this picture, 
there need not be a causal connection between star for- 
mation events in t he plane perpend icular to the large- 
scale flows (see also lElmegreenl [2OOOI ) . Rapid star forma- 
tion is a necessary requirement for this scenario to work. 

Numerical models of flow-driven cloud formation (we 
give only an early and t he most recent num erical 
wor k of each group, n amely iHe nnebelle fc Peraultj 



and iHennebelle et al.l II0O8: K ovama fc Inutsukaf 



1998 



200C 



and ^Ino ue &: Inutsukal I2OO81: iH eitsch et al. '2005' and 



Heitsch et al.l 2008bt IVazauez-Se madeni et al. 2006 and 
Vazquez- Semadeni et al.ll2007t ) have identified the ther- 



mal and dynamical instabilities that are responsible 
for the rapid fragmentation of the nascent cloud (see 
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IHeitsch et a l."2008a' for an assessment of the roles of the 
physical processes). Despite these promising successes, 
many questions about the physics at play remain unan- 
swered, among one of the most pressing is the role of 
magnetic fields during the cloud formation process. 

The role of magnetic fields in the flow-driven cloud 
formation scenario has been largely envisaged as one 
of "guiding the flows" to assemble the clouds, whether 
in form of the Parker instability along galactic spi- 
ral arms ()Parkeij [l966, 1967; Mouschovias et al. 1974), 
in a generally turbulent interstellar m edium (ISM) 
(jPassot et al.lll995HHartmann et al.ll2001h. or during the 
sweep- up of gas in s piral shocks |Kim fc OstrikeH [20061 : 
iDobbs fc Pried [20081 ). Based on the models of Passot et 
al., Hartmann et al. suggested that the field orientation 
with respect to the fiows selects the locations of cloud for- 
mation, namely that clouds will only form if the fields are 
aligned with the flows. A perpendicular field will reduce 
the compression of the post-shock gas and thus will limit 
the s trong cooling and the thermal instability (TI, Fieldl 
[19651) necessary for the rapid flow fragmentation and the 
build-up of high-density co ntrasts (jBurkert fc HartmannI 
[2003 [Heitsch et al.l[2008bl[a[ l. 

Given sufficiently high strengths, fields aligned with 
the inflows can suppress the dynamical instabilities re- 
sponsible for the generation of turbulen ce, namely the 



non-linear thin shell insta bility (NTSI, 



for a magnetic version see IHeitsch et al 



Vishniad [l99l 



2007) and the 



Kelvi n-Helmholtz instability (KHI, e.g. IChandrasekharl 
196ll more rece iitly[K unz 2008, and for numerical studies 
Palotti et al.ll2008l ). Yet magnetic fields arc intrinsically 



th ree-dimensional, and alr eady two-dimensional models 
bv llnoue fc Inutsukal (|2008l ) show that even for fields per- 
pendicular to the infiow, cold (albeit diffuse) clouds can 
form. 

Thus, three-dimensional models of flow-driven 
cloud formation i nclud ing magnetic fields are needed. 
IHennebelle et al.l (|2008f ) present a first approach to the 



problem, modeling the formation of a cloud in converg- 
ing, perturbed fiows, including fields and self-gravity. 
Here, we focus on the effects of magnetic field strength 
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and orientation on the early stages of flow-driven cloud 
formation. We work in the ideal MHD limit (i.e. we 
do not explicitly consider ambipolar drift or resistivity), 
and we do not include gravity in the models. 

All our models start out with field strengths below 
equipartition with the kinetic energy of the inflows. At 
a factor of 4.3 below equipartition - corresponding to an 
absolute field strength of 5^G at flow densities and ve- 
locities of 1 cm"'^ and 16 km s~^ - , the fields already 
suppress the dynamical instabilities (and thus the gener- 
ation of turbulence) leading to slab- like molecular clouds, 
while weaker fields - at 2.5/iG, corresponding to a factor 
of 17 below equipartition - lead to clouds closely resem- 
bling the hydrodynamical case, albeit with more coherent 
filaments. Fields at 0.5/iG perpendicular to the inflows 
suppress the build-up of massive clouds in the collision 
plane, while they lead to the formation of diffuse, cold 
filaments perpendicular to the mean field, remini s cent o f 
the cold m clouds discussed by Hcilc s fc Trolandl (|2003f ). 
A tangled field allows the assembly of substantial column 
densities in regions where the lateral field component is 
small or vanishing. Our results are consistent with the 
notion tha t magnetic fields selec t the location of cloud 
formation (jHartmann "eTalllMlh . 

2. TECHNICAL DETAILS & PARAMETERS 
2.1. Athena 

Calculations were performed with Athena 
([Gardiner & Stone 2005, 2008), an unsplit, second- 
order accurate Godunov scheme, using the c orner 
transport upwind metho d ([Gardiner fc Stonj|2008f ) and 
a linearized Roe solver (|Roelll98l[ ). The divergence of 
the magnetic field is kept zero by using constrained 
transport (Evans & Hawlej] 119881 ). Dissipative terms 
(viscosity, heat conduction and resistivity) are not 
explicitly included. For a det ailed description and test 
resul ts, th e reader is re f erred to lGardiner fc Stond ()2005l . 
[2008h and lStone etHl ([200l . 

We implemented heating and cooling as an additional 
energy source term at 2nd order in time. A tabulated 
cooling function provides the energy change rate as func- 
tion of density and temperature at each grid cell. We de- 
cided to keep the iterative approach we had used in our 
earlier studi es of cloud formation (jHeitsch et al.l l2005l 
l2006l l2008by with a slight modification. Instead of ad- 
vancing the fluid evolution at the usual time step given by 
the Courant-Friedrichs-Levy (CFL) condition and subcy- 
cling on the energy equation in case the cooling timescale 
is shorter than the CFL timestep, we lower the CFL 
timcstep according to 

At = AtcFL niin(l, {r./AtcFLT), (1) 

with < p < 1. For increasing p, small cooling timesteps 
will control the overall CFL timestep. The earlier ver- 
sion (see references above) would be equivalent to p = 0. 
Yet this choice can lead to inconsistencies in the hy- 
drodynamical evolution once the cooling timesteps get 
substantially shorter than the fluid timesteps, leading 
to a numerical overemphasis of the acoustic mode of 
the TI, since regions can cool substantially without ac- 
counting for the resulting pressure drop in the dynamics. 
While these inconsistencies may not affect the overall 
results, they turn out to affect the stability of the so- 
lution. For the models presented here, 0.5 < p < 0.7 



yields a stable and accurate solution. While the itera- 
tive approach is more time-co nsuming than a direct in- 
tegration as suggested by e.g. IVazguez-Semadeni et al.l 
(2007), it allows us to take into account the increas- 
ing thermal timescale of a fluid parcel cooling down to- 
wards the thermal equilibrium line (see e.g. Fig 3 of 
iHeitsch et al.|[2008al ). This effectively keeps a fraction of 
40-50% of the gas mass in the thermally unstable warm 
neutr al medium (WNM), con sistent with observed frac- 
tions (jHeiles fc Trolandll2003l: also gSS- 



2.2. Setup and Parameters 

The initial conditions and flow parameter are similar 
to those used in our earlier work (jHeitsch et al.ll2008b) . 
Two identical uniform flows collide head-on at an inter- 
face whose position is perturbed using an amplitude that 
is chosen randomly, but with the constraint that only the 
long wavelengths (down to a quarter of the box size) are 
non-zero. The maximum perturbation amplitude mea- 
sures 10% of the box size. The flows are otherwise in 
thermal equilibrium and do not have substructure. Thus 
we can test the most unfavorable conditions for substruc- 
ture and turbulence generation in the resulting clouds. 

The simulation domain has a resolution of 256^, and 
the linear box size measures 44 pc in all cases, result- 
ing in a nominal resolution of 0.18 pc. The boundaries 
in y (and z) are periodic, while those in x are set to 
a constant inflow speed of vq = 16 km s~^, a density 
of no = 1 cm^'^, a temperature of Tq = 8 x 10^ K 
and a magnetic field (if applicable). The adiabatic ex- 
ponent 7 = 5/3. The flow speed vq is consistent with 
shock speeds along streamlines p assing through G alactic 
spiral shocks at the solar circle (jShu et al.lll972[ ). Our 
thermal pressure is a factor 2 — 3 high er than the pres- 
sures estimated for the CNM /WNM () Jenkins fc TrippI 
[2OOII: iHeiles fc TrolM^dl [2001 . increasing the ratio of 
ther mal over magnetic press ure f3th above observed val- 
ues ()Heiles fc Trolandl l2005l hereafter HTG5). Yet the 
measure relevant for the dynamical importance of mag- 
netic fields is the ratio of the ram pressure over the 
magnetic pressure, (3ram, since the instabilities respon- 
sible for turbulence generation and fragmentation de- 
pend on the ratio of the flow velocities over the Alfven 
speed (Cliandrasckhar 1961; Hcitsch et al. 2007). In the 
diffuse, local ISM, kinetic and magnetic energy are in 
approximate equipartition {(3ram ~ 1; HT05), while in 
regions of ordered, large-scale flows - such as Galac- 
tic spiral shocks or expanding supernova shells - with 
uo ^ 10 km s~^, Pram > 1 should be expected. 

There are four classes of models: hydrodynamical (se- 
ries H), field aligned with flow (series X), field perpen- 
dicular to the flow (series Y), and (series XR) a uniform 
field component aligned with the flow plus a random fleld 
component of similar size, consistent with (although a 
little smaller than) observ e d magnetic field strengt h es- 
timates (e.g. iHeileal [19961: iBeckl l2004bl: Ha^lMffl . In 
the latter series, we do not perturb the collision interface 
but rely on the tangled field component to trigger the 
dynamical instabilities. Table [T] summarizes the model 
parameters. Self-gravity is not included in the models. 

To initialize the random field component, we set the 
amplitudes and phases of e.g. the x-component of the 
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TABLE 1 

Model Parameters 



Name 




ByO [mG] 




Pth 


0ram 


H 


0.0 


0.0 


0.0 


oo 


OO 


X25 


2.5 


0.0 


0.0 


4.5 


17 


X50 


5.0 


0.0 


0.0 


1.1 


4.3 


Y05 


0.0 


0.5 


0.0 


110 


430 


XR25 


2.5 


0.0 


2.5 


2.2 


8.5 



Note. — 1st column: model name, 2nd: magnetic field 
strength B^; 3rd: By, 4tli: random field Brms, Sth: thermal 
plasma f), 6th: ram plasma /3. 



(edge-centered) vector potential to 

max 

Ax{x,y,z)= ^ \k\'P sin{Kx + kyy + k:,z + (l)lj f.), 

(2) 

where \k\ = + ky + k'^ and e.g. kx = 2-Ki/Lx with the 
box length L^- We set p = 4, mimicking a (steep) turbu- 
lent energy spectrum as observed in detailed numerical 
simulations of magncto-hydrodynamic turbulence (e.g. 
[Cho & Lazarian 2003). The wavenumbers kx.y^z are cho- 
sen such that 1 < |/e| < 4, i.e. all combinations of (i, j, k) 
in the sum over fc-space are used that satisfy the con- 
straint on The phases (f)^^^^ i"^ A:-space are chosen 
from a uniform random distribution. Each vector po- 
tential component Ax y z requires a separate phase array 

This formula tion in real space instead of in Fourier 
space (see e.g. iMac Low et al.lllQQSl: iStone et alj 119981: 
iLemaster fc Stoni 120081 for velocity fields) allows us to 
easily regenerate the vector potential (and the field) at 
the inflow boundaries by 

Axi±Lx/2, y, z, t) = Ax{±{Lx/2 + VQt),y, z), (3) 

where the negative value refers to the lower x-boundary, 
and the positive to the upper one. The face-centered 
fields are then computed from the vector potential by 
B = V X A. 

The choice of the wave-number range 1 < < 4 does 
not constitute a restriction in terms of generality of our 
simulations, since the energy distribution over spatial 
scales is determined by the (steep) power law index p. 
This is fortunate in a sense, since the generation of the 
boundary conditions (eq. [3]) would consume substan- 
tially more time if wc had to sum over all available in 
equation 

2.3. Physical Interpretation of the Initial Conditions 

Obviously, our initial conditions are somewhat ideal- 
ized, e.g. generally, the flows would be expected to have 
substructure, the flows might not be expected to collide 
always head-on, and the magnetic fields will have par- 
allel and perpendicular components with respect to the 
inflows. Yet the initial conditions can be seen as idealized 
versions of different physical environments. 

The case of uniform fields aligned with the inflows 
(models X25, X50) could be identified with the sweep- 
up of material by an expanding supernova shell along an 
ordered background field, or with the collision of two ex- 
panding shells in such a field. The initial field strength of 
5/iG (model X50) is close to the lo cal median (tot al) field 
strength in the CNM (e.g. HT05: lTrola"^l2005f) . Using 



Nakano & Nakamura's l)1978f ) expression for the critical 

surface density Nc = B j^J A.-k^G above which gravita- 
tional collapse is possible under flux-freezing conditions, 
the swept-up clouds would reach approximately O.SiVc 
after 12 Myr, while model X25 {Bxq = 2.5/zG) would 
be marginally critical at the same time. We defer the 
discussion of the mass-to-flux ratio in the clouds to a 
subsequent paper including gravity. 

An ordered field aligned with the fiows plus a large- 
scale random component of similar amplitude (model 
XR25) introduces a large-scale shear and might be con- 
sidered a general situation for sweep-up of gas in spiral 
shocks, while the perpendicular field case (Y05) would 
address the (probably common) situation of an oblique 
field whose lateral component is amplified by flow com- 
pressions. 

We emphasize that while we attempt to address the 
extreme situations of field orientations, the finite size of 
our simulation domains cannot fully capture the effects of 
the magnetic field's boundary conditions. These will be 
set on larger scales than our local simulations can cover. 
In that sense, our results should be viewed as providing 
insight into magnetized cloud formation under idealized 
conditions rather than under physically realistic ones. 

2.4. A Comment on Resolution 

We decided to keep the resolution of our models con- 
stant, foregoing a resolution study in favor of a param- 
eter study. Resolut i on eff ects have been discussed by 
iHennebelle &: AuditI (|2007D . In addition, we have per- 
formed a systematic resolution study for two-dimensional 
cloud formation models (unpublish ed - the mod e ls are 
simil ar to the ones discussed by IHeitsch et al.l (|2005l 
l2006f )). covering a factor of 32 in spatial resolution (from 
256^ to 8192^ cells). As has been pointed out, the critical 
length scale to resolve is the cooling length of the ther- 
mal instability. If not resolved, the thermal instability 
will be partially suppressed. At parameters of the WNM, 
the cooling length is on the order of a parsec, while for 
the cold neutral medium (CNM), it drops to a fraction 
of a parsec. Thus, while more substructures should form 
with increasing resolution, we expect our models to fol- 
low the general evolution of the thermal and dynamical 
instabilities sufficiently accurately for our purposes. 

3. RESULTS 

3.1. Morphologies 

Figure [1] summarizes the morphological effects of mag- 
netic fields during the build-up of a cloud. From top to 
bottom, it shows logarithmic column density maps of the 
hydrodynamical model H, and the four MHD models X25 
through XR25. The three columns stand for projections 
along each coordinate axis, namely along the inflow (a;- 
axis, left), and perpendicular to the inflow (along y and 
z-axes, center and right). The maps of the MHD-models 
show polarization vectors which have been determined 
by integrating the density-weighted Stokes Q a nd U pa- 
rameters along the respe ctive line-of-sight (see IZweibell 
[Tool IHeitsch et al.l[200l . 

3.1.1. Field parallel to inflow 

Model H shows the strong fragmentation due to ther- 
mal and dynamical instabilities similar to the models dis- 




Fig. 1. — Logarithmic column density projections (in cm^^) along the three grid axes (left: along inflow, center and right: perpendicular 
to the inflow) for models H, X25, X50, Y05 and XR25 as indicated, at t = 9.5 Myr after flow collision. The mean field direction is denoted 
by the symbols in the panel labels {— ^,t,x). 
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cussed by iHeitsch et alj ()2008bD . Specifically, the large- 
scale initial perturbation triggers the NTSI, due to whose 
rapid growth some of the dense material has already 
reached the inflow boundaries. Viewed along the inflow 
(top left panel in Fig.[T]), the cold dense fragments appear 
clumpy rather than fllamentary. 

Introducing a magnetic fleld aligned with the flow 
(model X25, second row) suppresses fragmentation com- 
pared to model H. The face-on view {left) shows several 
large-scale coherent filaments with denser cores. The 
edge-on views (center and right) demonstrate that the 
magnetic field is not dynamically dominant. The polar- 
ization vectors are aligned with local structures. 

Increasing the magnetic field (model X50, third row 
from top) seems to suppress much of the fragmenta- 
tion. Specifically, the NTSI is only very weakly (if at 
all) present, since the magnetic field is strong enough to 
suppress the lateral momentum tra nsport necessary for 
triggering the NTSI (jHeitsch et al.| [2007). Nonetheless, 
the flows still fragment, albeit into a tight network of fil- 
aments (X50, left) instead of a few large, more clumpy 
and fuzzy structures (models H, X25). The suppression 
of the NTSI leads to the formation of a more or less 
coherent filament in the lateral projection (center and 
right column for model X50). Local shear modes lead to 
strong distortions of the field from its initial alignment 
with the inflow, as indicated by the polarization vectors 
which mostly trace out the mean background field. 

3.1.2. Field perpendicular to inflow 

The introduction of a field perpendicular to the inflow 
changes the morphology completely (4th row of Fig. [U 
model Y05), despite the by a factor of 10 weaker field (see 
Table [T]). The perpendicular field breaks the symmetry 
in the plane of the flow collision, leading to filaments 
perpendicular to the mean field direction (note that the 
mean field in the left panel of the 4th row of Fig. [1] is ori- 
ented horizontally). These filaments form due to motions 
along the field lines, but perpend icul ar to the incoming 
flows (see also IHeitsch et al.ll2007l and llnoue fc Inutsukal 
l2008l for two-dimensional models). The magnetic field 
suppresses one degree of freedom in the gas motions, also 
leading to lower column density contrasts than in models 
H and X50. The two lateral views of model Y05 exhibit 
another effect of the perpendicular field. Seen along the 
mean field direction, a large scale NTSI-drivcn mode is 
discernible, while the projection perpendicular to the in- 
flow and to the mean field (Y05 right) just shows a slab 
(albeit with substructure). In the former, the field lines 
are just shuffled around and contribute to the dynam- 
ics only via the pressure term in the Lorentz force, thus 
lowering the column densities and broadening the slab 
(compare to center panel of model X50). In the latter, 
the tension term of the Lorentz force prevents the growth 
of the NTSI. This is evidence for the presence of inter- 
change modes in t he NTSI, similar to e.g. the Rayleigh- 
Taylor instability (jStone fc Gardineij[2007[) . StiU, mate- 
rial is free to move along the field lines (and thus perpen- 
dicular to the inflows), leading to the formation of the 
filaments parallel to the inflows. 

The magnetic field perpendicular to the inflow resists 
compression, leading to a suppression of the thermal in- 
stability, which is also mirrored in the total mass budget 
of all models (Fig. [2|). Model X50 has the highest frac- 
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Fig. 2. — Total masses against time, below and above T = 100 K. 
The perpendicular field reduces the compression and thus lowers 
the cold gas mass, while the field aligned with the flow leads to 
higher compressibility. 

tion of cold gas, due to the strong guide field leading to 
a strong compression of the gas, while model Y05 shows 
the smallest cold mass fraction, because the lateral field 
resists compression by the fiows, and thus reduces the 
cooling rates. 

3.1.3. Tangled field 

The bottom row of Figure [T] shows the maps for model 
XR25, which starts out with a uniform field aligned with 
the fiow at 2.5yLtG and a random field component of equal 
magnitude. Although the (varying) lateral field compo- 
nents contain 10 times as much energy as the perpen- 
dicular field in model Y05, the fields do not suppress 
the formation of clouds with column densities in excess 
of 10^^ cm~^; a tangled field is substantially less effi- 
cient in preventing compression. Since there are regions 
where the field will be aligned with the flow, it leads 
to a selection effect in the sense that the clouds form 
at positions where the lateral random components of 
the fields are weakest over time and/or where bends in 
the fields determine th e position of cloud formation (see 
iHartmann et al.ll2001[ ). The resulting clouds are more 
isolated, with larger voids between them (bottom left 
panel of Fig. [1]). The side view (bottom center and right) 
exhibits a diffuse halo of thermally unstable gas, material 
which is caught in the tangled field between the bounding 
shocks and the dense cold gas. 

3.2. Dynamics 

The cold mass fractions of models H and X25 (Fig. ^ 
are slightly lower than that of X50, indicating that the 
developing turbulence due to the flow fragmentation is 
also broadening the slab. While this notion is already 
suggested by Figure [1] it is confirmed by comparing the 
rms velocity dispersion in the cold (T < 300 K) gas 
(Fig. [3]), and it can also be gleaned from a more detailed 
look at the laterally averaged pressure profiles (Fig. [3]). 

Shown is a time sequence of the pressure profiles along 
the X-axis (i.e. along the inflows) for all models. In the 
absence of gravity, the slabs are all overpressured by the 
ram pressure of the colliding flows (solid lines). At early 
times, all five models show a drop in kinetic pressure and 
an increase in thermal pressure in the collision region. 
The flows have not fully fragmented yet, and the cooling 
is not in full strength yet because of the still low densities. 
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Fig. 3. — Density- weighted rms velocity dispersion against time 
for all models. Fields parallel to the inflows seem to suppress tur- 
bulence in the cold gas. 

With evolving time, the thermal pressure peak for model 
H drops due to increasing cooling. 

This is markedly different for model X50, where the 
thermal pressure continues to be enhanced by a factor of 
more than 2 above that of the inflow. Also, the kinetic 
pressure drops, reaching an approximate equipartition 
with the thermal pressure. This is due to the strong mag- 
netic guide field, which "splits" the cloud into a network 
of dense filaments with low-density, high-temperature 
voids in between (see top panel of 2nd column of Fig. [Tl 
model X50). Figure [5] offers a different view of the same 
phenomenon, showing the pressure-density distributions 
for all four models, at t = 9.5 Myr. 

The high-temperature voids of model X50 show up at 
logn « 0.5 and logP « 4.3, while the high-density fil- 
aments sit all on the stable low-temperature branch of 
the thermal equilibrium curve at T w 40 K. Note that a 
substantial amount of the gas mass is actually thermally 
over-pressured, in contrast to model H. 

Reducing the field aligned with the flow (model X25) 
leads to pressure profiles (Fig. |4|) and thermal states 
(Fig. [5]) similar to the hydrodynamical model H. In other 
words, while the field in model X25 is non-negligible in 
the sense that its presence still makes a morphological 
difference (see Fig. [1]), it does not noticeably affect the 
overall dynamics of the cloud. 

The field is obviously dynamically important in model 
Y05. Because of the strong flow compression perpendicu- 
lar to the field lines, the magnetic pressure takes over the 
role of the thermal pressure, which leads to a substantial 
amount of thermally underpressured gas in model Y05 
(Fig. [4] and [5]). There is only a small amount of material 
at high (logn > 2) densities. 

Introducing the tangled field component on top of 
a field aligned with the inflows leads to an over- 
pressurization of the slab by a factor of more than 2 
(Fig. [4] right, model XR25). This is mainly due to a 
combined increase in magnetic and kinetic pressure, i.e. 
the tangled field leads to more turbulence than all other 
field geometries. The increase in kinetic pressure cannot 
be solely due to enhanced densities - model X50 should 
show a similar increase then. Although the tangled field 
in the diffuse gas phase is not force-free, it does not con- 
tribute perceptibly to turbulent motions in the inflows. 



as can be seen by comparing the kinetic pressure levels 
in the inflows between models XR25 and e.g. X50. Also, 
the kinetic pressure of model XR25 does not increase 
when moving closer towards the midplane, until one en- 
ters the post-shock region. The thermal pressure peaks 
in the diffuse envelopes due to warm gas being unable 
to cool down (see bottom right panel of Fig. [5]), but it 
drops back to the ambient thermal pressure at the cloud 
midplane (a; = 0). Obviously, the averaged thermal pres- 
sure alone is not a very accurate indicator of the cloud's 
physical state. 

Figure [6] summarizes the pressure evolution over time 
for all models. The pressures are sorted according to 
temperature regimes (rows; total, stable WNM, unsta- 
ble WNM, and CNM) and vector components (columns; 
total, X parallel to inflow, y perpendicular to inflow, 
and z again perpendicular to inflow). The total pres- 
sures {top row) stay constant with time except for the 
magnetic components of models XR25 and Y05 due 
to the compression of the transverse field component 
(see y column), and the y-component of the kinetic en- 
ergy of model XR25. Not surprisingly, the stable warm 
neutral phase does not show any strong pressure vari- 
ations with time. In the thermally unstable regime, 
300 < T < 3000 K, the models with transverse field 
components (XR25, Y05) show an increase in the y and 
z component of the magnetic pressure. A large fraction 
of the magnetic energy in model XR25 is thus stored in 
the diffuse cloud halo (see Fig. [T]). The magnetic fields 
also lead to a "redistribution" of kinetic energy between 
the vector components: the bulk of the kinetic energy 
of the cold gas in model H is found in the x-component, 
parallel to the inflow, indicating that the dynamics are 
dominated by the NTSI. The magnetic field models have 
at least as much kinetic energy in the transverse {y and 
z) components. For model XR25, the transverse com- 
ponents even dominate. The substantial overpressure at 
the midplane of model XR25 (see Fig. [5]) arises mainly 
from the transverse components of the magnetic and ki- 
netic energy in the cold gas. While for T < 300 K, the 
magnetic energies are larger than the kinetic ones, the ki- 
netic energy content dominates the total energy budget 
(top left panel of Fig. [6]). 

Figure [7] offers a simpler view of the field amplifica- 
tion depending on mean field direction. It shows the 
ordered (mean) field component and the random (tur- 
bulent) component against the time. Only the random 
components are growing with time, indicating that the 
fields are mainly amplified by fieldline stretching. Only 
models Y05 and XR25 have field components initially 
perpendicular to the infiows; these are growing linearly 
with time due to the compression by the flows. The ran- 
dom (or turbulent) fleld components of models X25 and 
X50 evolve in a similar way, indicating that a similar 
amount of energy is stored in the random component. 
For models X25 and Y05, the ordered and random com- 
ponents reach similar levels (see §4.3p . 



3.3. Turbulence and Thermal States 

The ratio of gas mass in the thermally unstable regime 
to the mass in the WNM can be compare d to observa- 
tional constraints. Heil es fc Trolandl (2003) flnd a frac- 
tion of 48%. Figure [5] shows the mass fractions for all 
models against the velocity dispersion in the CNM, av- 
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Fig. 4. — Pressure profiles along the inflow direction (x), averaged over the perpendicular directions (y, z), for three times as indicated, 
and for all models. Shown are total (solid line), kinetic (dashed), internal (dash-dot) and — if applicable - magnetic (dash-3-dot) pressures. 
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Fig. 5. — Greyscale-coded mass fraction of thermal pressure against volume density for all models as indicated in the plots, at t = 9.5 Myr. 
The solid line indicates the thermal equilibrium curve, while the vertical and horizontal dashed lines denote the initial conditions in density 
and pressure. Diagonal dashed lines stand for isotherms at T = 10*, 300 and 10 K as indicated. 



eraged over a time interval 7.5 < t < 9.5 Myr. Higher 
fractions of thermally unstable gas are found at higher 
velocity dis persions, a correlation al ready p ointed out 
earlier (e.g. iGazol et al.ll2001l:IAudit fc Hcnne bellell2005l : 
iHeitsch et al.l l2006t iHennebelle fc Audit. i2007f) . How- 
ever, the magnetic field introduces a second dependency: 
of the three models with the highest velocity dispersion 
(H, Y05 and XR25), model Y05 has the highest frac- 
tion of thermally unstable gas. As discussed above, in 
this model the lateral field components reduce the den- 
sity contrasts, keeping a substantial fraction of the gas 
in the thermally unstable regime. The large velocity dis- 



persion of model Y05 arises from the large-scale NTSI 
interchange mode (see Fig.[T]). 

It is maybe not surprising that the structures which 
bear closest resemblance to t he CNM clouds studied by 
iHeiles fc Trolandl (|2003ll2005[ ) belong to the model with a 
thermally unstable mass fraction closest to the observed 
one (Y05). 

4. DISCUSSION 

4.1. The Role of Magnetic Fields for Cloud Formation 

The field strength, and the orientation of the mean 
magnetic field with respect to the flows sweeping up the 
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Fig. 8. — Mass fraction of thermally unstable gas (300 < T < 
3000 K) with respect to all WNM (T > 300 K) against the velocity 
dispersion in the CNM (Fig.|3]l, averaged between 7.5 and 9.5 Myr. 
The error bars show errors on the mean. 



gas play a crucial role for the flow-driven formation of 
molecular clouds (Fig. [T]). If the fields are dynamically 
dominant, the only chance to build up substantial clouds 
is by channeling the flows along the fields. This is the 



situation shown in model X50, and it is also borne out by 
simulation s of molecular c loud formation in galactic spi- 
ral arms (Kim fc Ostrikeij[2006.) . where the clouds tend 
to be oriented perpendicularly to the large-scale field 
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(also possibly visible in the models by iDobbs fc Pried 
|2008| ). until sufficient material has been accumulated 
that they decouple dynamically from the large-scale field. 
Similarly, for the sweep-up of material by e.g. Hll-regions 
or supernova shells, one would expect the densest clouds 
to appear at the locations where the field is perpendicular 
to the shell (see Fig. 2 of iHanavama fc Tomisakal l2006l . 
although the effect might be less cle ar in a highly tu rbu- 
le nt environment, s ee Fig . 10 of Bals ara et al.|[200^ . 

iHartmann et al.l (|2001D po int out - based on simula- 
tions bv iPassot et al.l (|1995( ) - that dynamically weak 
(but not necessarily ordered) fields would lead to a gen- 
eral selection effect for the formation of molecular clouds. 
Since /3ram > 1, the flows stretch out the fieldlines, lead- 
ing to a natural alignment. In this picture, clouds form 
in the bends of large -scale fields (see Figs. 4 & 5 of 
IHartmann et al.ll200lD . Such a scenario is to some ex- 
tent addressed by model XR25, where varying field ori- 
entations entail a local selection effect, picking out the 
formation sites of molecular clouds over e.g. a broad 
shock front. Note that while the field is dynamically 
weak (Pram = 8.5, /3th — 4.3) in the initial conditions 
(and in the infiows) of model XR25, Pth < 1 within the 
cloud (Fig. SD. 

This selection effect comes about because already a 
small oblique component can be amplified sufficiently to 
withstand the compression, preventing the high de nsities 
needed for cloud formation (|Hartmann et al.ll2001| ) . Such 
a situation is addressed in the extreme by model Y05. A 
field perpendicular to the sweeping-up fiow can suppress 
the formation of massive clouds, although the three- 
dimensional situation is much less clear-cut than its one- 
dimensional counterpart (see e.g. iMcKee fc HollenbachI 
\Tmi iBergin et al.l[20nl . In one dimension, a density 
increase by a factor of 100 from e.g. n = 1 cm~^ to 
100 cm~'^ would entail the same factor for the magnetic 
field strength since B (x n. Figure [9] shows this is not 
the case in three dimensions. The weak perpendicular 
field (model Y05) has been amplified by a peak factor 
of « 30, while the density has increased by up to a 
factor of 300. Generally, our models show a weak cor- 
relation of field strength with density over the whole 
thermal range, from the WNM to the CNM, consis- 
tent with observati ons of the field-density r elation in the 
WNM and CNM (iTroland fc H eiles' T98 l; HT05), and 
with numerical results ( e.g. | d~A villcz fc Breitschwerdtl 
120051 iHennebelle et all l2008l ). For models X50 and 
X25 a weak correlation between field and density is 
not overly surprising. For model Y05, the decorre- 
lation^ is a consequence of the fact that material is 
still free to move alo ng the field lines p erpendicular 
to the original inflow (Hei tsch et al.l [2b07[ ). thus lead- 
ing to the build-up of higher-density fllaments per- 
pendicular to the fleld (but aligned with the inflow), 
see Figure [TJ Also, other effects, such as the ac- 
celeration ofniagne^ by turbulence 
(iZweibell 120021: iF'atuzzo fc A dams' 2002; 'Heitsch ct al.' 
120041 for ion-neutral drift, [Lazarian fc Vishniac 1999 for 
reconnection), or a decorrelation due to MHD waves 

^ The seemingly strongly correlated B(n) for n < 1 cm""' in 
model Y05 does not affect the argument. These are a few regions 
(low mass fraction) at the edges of the expanding slab, subjected 
to numerical reconnection. 



(jPassot fc Vazguez-Semadenil I2003D could explain the 
observed weak correlation. 

Magnetic fields will rarely be completely uniform. 
Model XR25 tests the more general case of a uni- 
form field at 2.5/iG and a random component of equal 
size, consistent with (although slightly lower than) ob- 
servational estir nates for m a gnetic field strengths i n 
the diffuse gas fHeilef^ 1^91 [Bec^ l2004bl: lHaJ[2006l) . 
iGiacalone fc Jo kipii (20()3) showed in a two-dimensional 
numerical experiment that pre-existing perturbations in 
the inflows can lead to substantial magnetic field ampli- 
fication due to a rippling of the shockfront and subse- 
quent fieldline stretching. We observe a similar effect in 
model XR25, although our Mach numbers are substan- 
tially lower (their study address ed the pr opagation of a 
supernova shock front). Henn ebelle et al.l (2008) perturb 
the velocities of the inflows and find only a modest in- 
crease of the field strength. Clearly, the initially tangled 
field leads to rather different dynamics in the forming 
cloud (FigsIHEl). 

Models X25 and X50 demonstrate that not only the 
field orientation will play a role during cloud formation 
(see model Y05), but also the field strength, since all 
the instabilities involved have threshold limits for the 
field strength - at least in two dimensions. It might well 
be that the stronger field in model X50 suppressing the 
formation of filaments could be offset by higher inflow 
speeds or substructures in the flows. This remains to be 
studied. 

4.2. Turbulence and Thermal States 

Fields aligned with the inflows tend to suppress the 
NTSI, and thus lead to an approximate equipartition 
between the spatial components of the kinetic energy 
in the cold and in the thermally unstable gas (Fig. Hj 
bottom couple of rows). For comparison, the hydro- 
dynamical model H has the bulk of the kinetic energy 
in the (flow-aligned) x-component. Magnetic fields may 
play an important role to isotropize highly directional 
flows. Thus, searching for observational signatures of 
flow-driven cloud formation should focus on the warm, 
diffuse gas phase, since the inflow signature will be erased 
in the cold dense gas. 

Fractions of thermally unstable gas (Fig. [5]) for flow- 
aligned fields (models X50, X25) are lower tha n val- 
ues o bserved for diffuse CNM clouds (Heilcs fc Tr olandl 
|2003|) . A lateral field component results in a thermally 
unstable gas fraction of « 50%, consistent with observa- 
tions. Based on these findings, one could feel tempted 
to extend the above argument about the selection ef- 
fect introduced by magnetic fields: not only could mag- 
netic fields control the locations of molecular cloud for- 
mation, but they also could lead to "failed" molecular 
clouds, i.e. diffuse atomic hydrogen clouds, if there is 
a non-negligible field co mponent perpen dicular to the 
sweeping-up fiow (see also llnoue fc Inutsu ka (2008) for a 
similar argument based on two-dimensional simulations) . 

4.3. Ordered vs. Random Component 

Another observational constraint is given by the ra- 
tio of the ordered over the unordered (or turbulent) field 
component. The observational evidence poi nts to the 
components being of simi l ar magnitu de (e.g. ,Heileslll996l : 
lBeckll2004bl: lBeckll2004al: iHanl 120061: see also discussion 
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Fig. 9. — Magnetic field strength against volume density for models as indicated in the plots. The colors denote temperatures, and the 
intensity the mass fraction. Generally, there is no clear correlation between field strength and density. The steeper of the two dashed lines 
denotes B oc n, the flatter one B oc n^/^. Dotted lines denote the initial conditions. 



in HT05). A direct comparison to our models is ham- 
pered by the fact that in order to see the varying com- 
ponent, sufficiently large scales need to be addressed, 
which is why HT05 argue that their observed median 
field strength of 6/iG should be identified with the to- 
tal magnetic field strength. Likewise, it is not obvious 
that the components should be of equal magnitude lo- 
cally everywhere. Bearing this limitation in mind, it is 
clear from Figure [7] that only for models X25 and Y05 
the components are comparable. 

4.4. Gravity 

We deliberately left out self-gravity in our simulations, 
in order to get a clearer view of the role of magnetic 
fields during the early cloud formation phase. Thus, our 
clouds are only confined by the ram pressure of the in- 
flows, and at later stages, the dynamics of the clouds are 
probably underestimate d since gr avity as a source of tur- 
bulence is missing (e.g. iField et al. .2008,) . As a result of 
the restricted physics, a comparison of our models with 
observations is only meaningful for models where gravity 
is not expected to play a role, i.e. for model Y05 ad- 
dressing the formation of diffuse HI clouds. For all other 
models, we expect gravity to be relevant during the cloud 
formation process (jHeitsch fc Hartmannll2008[ ). 

5. SUMMARY 



Extending our previous work and complemen ting a 
mode l by Hennebelle et al. (2008; see also Baneriee et al.l 
l2008t ). we study the role of magnetic field strength and 
orientation on the process of flow-driven cloud formation. 
Our models include the usual heating and cooling effects, 
allowing rapid fragmentation of the flows, they use uni- 
form inflows to study the most unfavorable conditions 
for structure formation, and they envisage the formation 
of clouds in two head-on colliding flows, i.e. the extreme 
case for building up massive clouds. We do not include 
self-gravity, focusing on the early stages of cloud forma- 
tion, during which gravity might be less important. 

Under these conditions, we find that the effects of mag- 
netic fields on the morphology and on the thermal state 
of the resulting clouds depend very strongly not only 
on the field orientation with respect to the inflow, but 
also on the fleld strength. Initial field energies are be- 
low equipartition with the kinetic energies (by a factor 
of 4.3, corresponding to a field strength of 5fj.G for our 
flow parameters) even for the strongest field case in our 
study (model X50), yet they result in significantly differ- 
ent cloud properties than those for a field weaker by a 
factor of 2 (model X25, 2.5/iG). Magnetic fields also lead 
to a redistribution of the inflow energy to the transverse 
spatial directions (Fig.[6|). Hence searching for signatures 
of colliding flows should focus on the diffuse gas phase, 
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since the cold gas will have no memory of the original 
flow direction. 

Not surprisingly, weak magnetic fields (0.5/iG) perpen- 
dicular to the inflows can suppress the build-up of mas- 
sive clouds (model Y05). Yet substructure still can arise 
in the post shock gas, in the form of diffuse filaments per- 
pendicular to the field, and of wave-like patterns (pos- 
sibly magnetosonic waves). The filament s are a conse- 
quence of lateral gas transp ort (see also iHeitsch et al.l 
I2007t llnoue fc Inutsukal '20081 . The straight-forward cor- 
relation B ot: n is not obeyed (Fig. [9]). Mass fractions of 
thermally unstable gas for the model with a lateral field 
component (Y05) are consistent w ith observed values for 
diflfusc HI clouds (|Heiles fc Trola nd 2003). For aU other 
models, the fractions are lower (Fig. [8]). The ratio of 
ordered vs. random field component is consistent with 
observations only for the weak-field model X25, and for 
the diffuse HI cloud model Y05 (Fig. [3). 

A weak (2.5/iG) uniform field together with a ran- 
dom component of equal size leads to a strong over- 
pressurization of the cloud due to a combined increase 



of magnetic and kinetic pressure (Fig. [4]) , with the mag- 
netic pressure dominating the thermal pressure within 
the cloud. High column densities are assembled at loca- 
tions where the perpendicular field component is weakest 
over time. Thus, a tangled field can lead to a selection 
effect for cloud formation while not preventing it globally. 

Our numerical models address the ideal MHD limit, 
i.e. we do not take into account ion-neutral decoupling or 
resistive dissipation. It remains to be seen how non-ideal 
MHD processes affect the struct ure formation during the 
build-up of the clouds (e.g. Jnoue fc InutsukalfeOOSl . 
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